Elliptic Flow and Dissipation in Heavy-Ion Collisions at E l8Lh ~ (1—160) A GeV* 
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Elliptic flow in heavy-ion collisions at incident energies i?i a b — (1—160)^4 GeV is analyzed within 
the model of 3-fluid dynamics (3FD). We show that a simple correction factor, taking into account 
dissipative affects, allows us to adjust the 3FD results to experimental data. This single-parameter 
fit results in a good reproduction of the elliptic flow as a function of the incident energy, central- 
ity of the collision and rapidity. The experimental scaling of pion eccentricity-scaled elliptic flow 
versus charged-hadron-multiplicity density per unit transverse area turns out to be also reasonably 
described. Proceeding from values of the Knudsen number, deduced from this fit, we estimate the 
upper limit the shear viscosity-to-entropy ratio as ?7/s ~ 1 — 2 at the SPS incident energies. This 
value is of the order of minimal rj/s observed in water and liquid nitrogen. 
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I. INTRODUCTION 

The elliptic flow (v?) of produced particles is one of 
the most sensitive obscrvablcs which brings information 
about the degree of collectivity during the expansion 
stage of heavy- ion collisions [l], 0, H[ . When the collec- 
tivity is strong, like in the case of ideal hydrodynamics, 
the elliptic flow takes the highest value (the so called 
hydrodynamic limit). If the collectivity is weak, like 
in a dilute system of weakly interacting particles, it is 
close to zero. Therefore, it is not surprising that v% 
is highly sensitive to dissipative effects which are as- 
sociated with attenuation of the collectivity during the 
expansion stage. The elliptic flow turns out to be con- 
siderably reduced by the cascade "afterburner" follow- 
ing the hydrodynamic freeze-out [j, [E 0] as well as by 
viscosity effects @, H i, EE OH UM- Note that the af- 
terburner can be considered as a strong viscosity effect 
at the final stage of the fireball expansion. 

All above mentioned analyses of the dissipative effects 
were carried out for high incident energies in the RHIC 
region. In this paper we would like to turn to lower 
energies of the 

SIS-AGS-SPS region. The experimental data [U, 

Hi HE HE mm QaliS 

in this energy region are much more fragmentary 
than in the RHIC domain. In the future they will 
be essentially complemented by new facilities, FAIR 
in Darmstadt and NICA in Dubna, as well as by 
experiments within the low-energy-scan program at 
RHIC. The available data were analyzed within kinetic 
M, l2l HI HI and 3FD M, HI models. It was found 



that the 3FD approach noticeably overestimates the 
data at the SPS energies independently of the stiff- 
ness of the equation of state (EoS) and stopping power 
used by the model. Recently the SPS data were studied 
within a hybrid hydrocascade model [27j which includes 
afterburner effects. It was demonstrated that the after- 
burner indeed essentially reduces the hydrodynamic t>2 
values which, however, still exceed their experimental 
values. 

In fact, there were no studies of dissipative effects re- 
lated to the elliptic flow in the AGS-SPS energy range. 
In this paper we would like to present such a s tudy 
based on simulations within the 3FD model [25Ll26Ll28j|. 
In order to estimate the dissipation, we use the ap- 
proach suggested in Refs. [IE [3(1 El]. We also study 
scaling of the elliptic flow with the midrapidity density 
of produced charged particles, which was experimen- 
tally revealed in Refs. (3II32I]. 



II. 3FD MODEL 

In Ref. p8l | we have introduced a 3-fluid dynamical 
model for simulating heavy-ion collisions in the cm. en- 
ergy range ^i^v = 2.3-30 GeV 1 (or £ lab = (l-500),4 
GeV in terms of lab. energy of the beam), which over- 
laps with the SIS-AGS-SPS energy range and covers 
the domains of future FAIR and NICA facilities. The 
3FD model is a straightforward extension of the 2-fluid 
model with radiation of direct pions [H, [H| and the 
(2+l)-fluid model (3f|. These models were extended in 
such a way that the created baryon-frec fluid (a "fire- 
ball" fluid) is treated on equal footing with the baryon- 



*This paper is dedicated to our long-term coauthor and friend, 
V.N. Russkikh, without whom this and many other works would 
be impossible, 
te-mail: Y.Ivanov@gsi.de 



The upper limit is associated with computational demands 
of the 3FD code. The lower limit cuts off the region where 
applicability of the hydrodynamics becomes questionable. 
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rich fluids. For the fireball fluid we have introduced a 
certain formation time, during which it evolves without 
interaction with other fluids. 

The input required by the model consists of the EoS 
and inter-fluid friction forces. Our goal is to find an 
EoS which is able to reproduce in the best way the 
largest body of available observables. The friction forces 
determine the stopping power of colliding nuclei and 
thereby the rate of thermalization of produced matter. 
In principle, the friction forces and the EoS are not in- 
dependent, because medium modifications, providing a 
nontrivial EoS, also modify elementary cross sections. 
However, currently we have at our disposal only a rough 
estimate of the inter-fluid friction forces [36[. In the 
present version of the 3FD model these forces are fit- 
ted to the stopping power observed in proton rapidity 
distributions. 

We have started our simulations [28[ with a simple, 
purely hadronic EoS [37| which involves only a density 
dependent mean-field providing saturation of cold nu- 
clear matter at normal nuclear density no = 0.15 fm 
with the proper binding energy —16 MeV and a certain 
incompessibility K. It includes 48 lowest mass hadronic 
states ,281] . This EoS is a natural reference point for any 
other more elaborated EoS. 

The 3FD model turns out to reasonably reproduce 
a great body of experimental data in a wide energy 
range from AGS to SPS. Figure [1] illustrates the re- 
sults of Ref . [28| . Here we use the compilatio n l38l o f 
experimental data from Refs. [H, E EE S EaE|. 
These data slightly differ in degree of centrality. For 
Au+Au collisions at AGS energies and Pb+Pb reac- 




FIG. 1: (Color online) Incident energy dependence of midra- 
pidity yield of various charged hadrons produced in central 
Au+Au and Pb+Pb collisions. The 3FD calculations are 
done with the hadronic EoS (K = 190 MeV). The compila- 
tion of experimental data is taken from Ref. [38l |. 



tions at E\ a b = 158 A GeV we perform our calculations 
taking the fixed impact parameter 6 = 2 fm. For Pb+Pb 
collisions at lower SPS energies we use b = 2.5 fm. The 
overall description of the data is quite good in the whole 
energy range under consideration. However, a certain 
underestimation of the kaon yield at lower SPS energies 
prevents the model from reproducing the " horn" in the 
K + /n + ratio observed experimentally in Refs. (45l. [46|. 

In the present simulations we use slightly different 
set of parameters as compared with that of Ref. [28| . 
This is caused by several reasons. We excluded all con- 
tributions of weak decays into hadronic yields, as it 
required by the NA49 experimental data. Strong de- 
cays of baryon resonances were updated. Scalar me- 
son /o(600) was included in the list of produced mesons. 
Some bugs in the code were corrected. All this inspired 
a refit of the model parameters. The best EoS, which 
reproduces the main body of experimental data, cor- 
responds now to K = 190 MeV (instead of K = 210 
MeV in Ref. [H). In fact, the EoS's with K = 190 
and 210 MeV are very close to each other. The en- 
hancement factor of the friction forces estimated from 
the proton-proton cross sections was reduced: we take 
the coefficie nt B h = 0.5 instead of (3h — 0.75 in Eq. 
(39) of Ref. [2g|. With this updated set of parameters 
the reproduction of available data is approximately the 
same as in [28j. 



III. ELLIPTIC FLOW 

The elliptic flow, defined as v-2 = (cos 2(f)) 0], is the 
second coefficient in Fourier expansion of the azimuthal- 
angle 2 dependence of the single-particle distribution 
function of a hadronic species a, (Ed 3 N a /d 3 p), 

j d 2 PT [(pi - pI)Ip\\ Ed 3 N a /d 3 p 

v ( 2 a) (y) = J - T , (1) 

/ d 2 p T Ed 3 N a /d 3 p 

where px is the transverse momentum of the particle, 
p x and p y are its x and y components. In calculating V2 
for pions and protons we take into account contributions 
of resonance decays. 

Figure [2] summarizes the 3FD results [H [H, [28[ for 
the elliptic flow. The calculations have been done for 
Au+Au collisions at b = 6 fm (SIS and AGS energies) 
and Pb+Pb collisions at b = 5.6 fm (SPS energies). It 
is seen that the proton data are well reproduced at low 
bombarding energies, where the squeeze-out effect dom- 
inates. The latter is caused by shadowing of expand- 
ing participant matter by spectator parts of the col- 
liding nuclei. Note that late freeze-out turns out to 



2 rfi is the angle with respect to the reaction plane. 
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FIG. 2: (Color online) Elliptic flow of protons (top panels) and pions (bottom panels) at midrapidity as a function of incident 
energy in mid-central Au+Au (at SIS and AGS energies) and Pb+Pb (at SPS energies) collisions. The 3FD calculations 
with soft (K = 130 MeV), intermediate (K = 190 MeV) and hard (K = 260 MeV) EoS's at standard freeze-out energy 
density £f rz = 0.4 GeV/fm 3 are displayed in right panels. Calculations with intermediate (K = 190 MeV) EoS and different 
freeze-out energy densities £f rz = 0.2, 0.4 and 0.6 GeV/fm 3 are presented in left panels. Compilation of experimental data 
is from Ref. fl5f | - The FOPI pion data are from [l^]: filled symbols correspond to positive pions, and open symbols, to 
negative pions; the rapidity average is taken over the interval —1.8 < y — y c .m. < 0. 



be preferable for pion data at low energies, while the 
reproduction of these data is still far from being per- 
fect. At higher energies, when the standard collective 
mechanism of the elliptic flow formation starts to work, 
the 3FD model noticeably overestimates both the pro- 
ton and pion data. The calculated results only weakly 
depend on the stiffness of the EoS. Tuning the freeze- 
out condition does not help very much to reduce the 
disagreement. The proton elliptic flow turns out to 
be quite insensitive to this condition. The sensitivity 
of the pion elliptic flow is higher. However, even very 
early freeze-out, which looks preferable for pions, does 
not allow to fit the data. The different sensitivity of 
the proton and pion V2 to the freeze-out condition is a 
consequence of three-fluid nature of our model. Satu- 
ration of i>2 in the fireball fluid occurs later than that 
in the baryon-rich fluids. Since contribution of the fire- 



ball fluid into the pion yield is larger than for protons, 
the saturation of the total pion i>2 happens also later. 
Variation of the inter-fluid friction does not improve 
agreement with the data either, while destroying the 
description of other observables. 

The calculations show that for e trz < 0.2 GeV/fm 3 , 
both the proton and pion elliptic flows stay practi- 
cally unchanged. Therefore, the results obtained for 
£frz = 0.2 GeV/fm 3 can be naturally associated with 
hydrodynamic limit 3 of the elliptic flow. 

It is natural to associate the overestimation of ex- 
perimental i>2 values with dissipative effects during the 



3 Note that these results correspond to the 3FD initial conditions. 
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FIG. 3: (Color online) Same as Fig. [2] but for 3FD re- 
sults reduced in accordance with Eqs. ([2}, Q with c 3 a tI — 
2.3 mb. The results of 3FD calculations for the intermedi- 
ate (K = 190 MeV) EoS and data are shown for central, 
midcentral and peripheral collisions. Different experimental 
points (NA49) for the same incident energy and centrality 
correspond to different experimental methods of vi determi- 
nation, see Ref. [l{J. 



expansion and freeze-out of the participant matter. To 
take these effects into account we use an empirical for- 
mula suggested in [3(| : 



hydro 



''2 



1 + Kn/Kno 



(2) 



Here V2 is the observed value of the elliptic flow, ii)/^ 10 
is its hydrodynamic limit, Kn is an effective Knudscn 
number defined as 



Kn 



A 
R 



(3) 



where A is a mean-free path of a particle and R is a 
characteristic size of the system (e.g., the radius of the 



nucleus), Kno ~ 1 is a constant. The recent transport 
calculation in two spatial dimensions (3lj resulted in 
Kno — 0.7. We use this value of Kno for our estimates 
below. 

As argued in Ref. [3(J, the Knudsen number can be 
represented in the form 



1 C s CT tr dN t , 



Kn 45 dy 



(4) 



Here dN tot / dy is the total (charged plus neutral) hadron 
multiplicity per unit rapidity, which equals approx- 
imately 3/2 of the charged-hadron rapidity density 
dNch/dy, cr tr is the transport cross section, c s is the 
sound velocity in the medium, and S is the transverse 
overlap area between two colliding nuclei. The latter is 
defined as S = ny / (x 2 ) (y 2 ) (49[, where (x 2 ) and (y 2 ) 
are mean values of x 2 and y 2 over the overlap zone 4 
defined by the collision geometry. 

Similarly to Ref. we fit the data displayed in 

Fig.Uby using Eq. © with Kn defined by © and v^ ydl '° 
taken from the 3FD calculation with the late freeze- 
out, i.e. at £f rz = 0.2 GeV/fm 3 . The rapidity density 
dN c h/dy is also calculated within the 3FD model. Note 
that the charged-hadron rapidity distributions are well 
reproduced by the 3FD model (see Fig.Q]). The reduc- 
tion factor (l + Kn/Kno) -1 is applied only at i?i a b > 4A 
GeV, since it is not justified at lower energies due to 
importance of the squeeze-out effects. Indeed, Eq. ([2]) 
was deduced from simulations of unbiased expansion 
of a system into the transverse directions (30j . In this 
case the elliptic flow can be only positive. The squeeze- 
out means that the transverse expansion is screened by 
spectators. The latter results in a suppressed and even 
negative elliptic flow. Therefore, at lower incident ener- 
gies purely hydrodynamical results for the late freeze- 
out (ef rz = 0.2 GeV/fm 3 ) are presented. Probably, the 
elliptic flow in the squeeze-out region should be cor- 
rected for the dissipativc effects but certainly not by 
means of Eq. ([2]) . 

The results of our fit are presented in Fig. [3] One can 
see that a good reproduction of the AGS and SPS data 
is achieved with the help of a single fitting parameter 
c s cr tr ~ 2.3 mb. For the midcentral collisions this leads 
to the estimate Kn ~ 0.7 at the midrapidity. Intro- 
duction of the dissipative correction also substantially 
improves the rapidity dependence of the elliptic flow. 
This is indeed seen from Fig. [3J Taking c 2 ~ 0.15 [5fj |. 
we arrive at the estimate Ct r ~ 6 mb for top AGS and 
SPS energies. This is a surprisingly low value of the 
cross section. However, taking into account that cr tr is 
the transport cross section (which is, in general, lower 
than the total one) and having in mind uncertainties in 



4 We follow the conventional definition of S [49 
from that in Refs. [30I . [3l{ by the factor of 4. 
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FIG. 4: (Color online) Elliptic flow of protons (upper panels) and charged pions (lower panels) in mid-central Pb+Pb 
collisions at = 158^4 GeV (left panels) and 40A GeV (right panels) as a function of rapidity. The 3FD calculations are 
performed at b = 5.6 fm with the intermediate EoS. Experimental data obtained by different methods are displayed: by 
the standard method («a(st)) and by the method of n-particle correlations (v2(n)). Full symbols correspond to measured 
data, while open symbols are those reflected with respect to the midrapidity. Circles in left panels show the updated NA49 
data [13] with the acceptance 0.05 < pr < 0.35 GeV/c for pions and 0.6 < pr < 2.0 Gev/c for protons. 



the Knudsen number definition this value seems to 
be still acceptable. 

Let us discuss now an app roximate scaling behavior 
of i'2 proposed in Refs. @,|32|. This behavior is observed 
when «2 scaled with the initial eccentricity e is plotted 
as a function of dN^/dy scaled with the cross section 
of the nuclear overlap S. Both e and S are determined 
by the collision geometry. In particular, e is defined as 

Below we calculate (x 2 ) and (y 2 ) with either the 
wounded-nuclcon (WN) or the binary-collision (BC) 
weights, for details sec Ref. [49]. These calculations arc 
based on the usual Woods-Saxon profile of nuclear den- 



sity. Within the 3FD model the initial nuclei are rep- 
resented by sharp-edged spheres. To be consistent with 
the model, we first calculated e and S using the Woods- 
Saxon parametrization with zero diffuseness (d = 0) and 
the BC weight. The obtained results are shown in Fig. [5] 
by the solid lines. 

Our scaling analysis is summarized in Figs. [5] and [7] 
First of all, we have found that there is no scaling for 
proton elliptic flow: points corresponding to different 
energies and impact parameters populate a relatively 
broad band rather than a universal line. This is not sur- 
prising, since scaling either with the BC or WN weights 
is inappropriate for nuclcons. Indeed, their rapidity 
spectra arc strongly constrained by the conservation of 
the baryon number. As a result, contrary to pions, the 
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b [fm] 

FIG. 5: (Color online) The transverse overlap area S and 
spatial eccentricity e as functions of impact parameter in 
Au+Au collisions for different surface diffusenesses of the Au 
nucleus (d) and different weights of averaging: the wounded- 
nucleon (WN) and the binary-collision (BC) weights [4j|. 



number of participating nuclcons is determined mainly 
by the initial geometry and it is not proportional to the 
number of binary collisions. 

At the same time the pionic «2 exhibits an approx- 
imate scaling behavior already for pure hydrodynamic 
calculation, see the top panel of Fig. [5] This scaling 
is more pronounced at higher densities of charged par- 
ticles, i.e. at higher incident energies. However, this 
scaling still substantially differs from the observed one, 
cf. experimental points in Fig. [HI When we correct 
the hydrodynamic results according to Eq. ([2]) (see 
the bottom panel of Fig. 0, the resulting V2 reveals 
a scaling behavior which is closer to the experimental 
results. By thin lines in Fig. [5] we also show the hy- 
drodynamic limits of the elliptic flow. The latter are 
obtained from the 3FD calculation with the late freeze- 
out (e frz = 0.2 GeV/fm 3 ). 

Although the scaling factors e and S, used in Fig. [5] 
(see the "BC: d = 0.0" curves in Fig. 0, are consistent 
with assumptions of the 3FD model, they differ from 
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FIG. 6: (Color online) Pion elliptic flow at midrapidity di- 
vided by the eccentricity as a function of the charged-hadron 
rapidity density in Au+Au (AGS energies) and Pb+Pb (SPS 
energies) collisions at different centralities. Bold lines in 
upper and bottom panels display, respectively, the hydro- 
dynamic vi calculated with the standard freeze-out energy 
density (£f rz = 0.4 GeV/fm ) and the hydrodynamic V2 cor- 
rected accordingly to Eq. (J5J). Thin lines represent the un- 
corrected hydrodynamic calculation with the late freeze-out 
(ef rz = 0.2 GeV/fm 3 ). The 3FD calculations have been done 
with the intermediate EoS (K = f90 MeV). Compilation of 
experimental data is from Ref. 



those used for scaling of the experimental data [3j . The 
experimental points in Figs. O and [7] are scaled with e 
and S calculated with the WN weight and d = 0.535 fm 
(49j (see the "WN" curves in Fig. |5J) . It turns out that 
e and S calculated with the BC weight and d = 0.5 fm 
agree fairly well with the "experimental" scaling factors. 

Upon applying the "experimental" scaling factors to 
the corrected elliptic flow, we obtain the scaling-like be- 
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FIG. 7: (Color online) The same as in bottom panel of Fig. 
[5] but scaled with "experimental" scaling quantities S and 
s. Elliptic flow at incident energies Eub > 10^4 GeV is 
displayed by bold lines, while that at lower energies < 
WA GeV - by thin lines. 



ticcably affect the elliptic flow in semicentral collisions, 
especially at high transverse momenta and marginal ra- 
pidities. However, for the midrapidity region in midcen- 
tral collisions, considered here, the effect of fluctuations 
is quite moderate. 



IV. DISSIPATION 

Having determined the effective Knudsen number, we 
can now estimate the role of dissipative effects during 
the expansion stage of a nuclear collision. Let us first 
attribute all this dissipation to the fluid viscosity. To es- 
timate the latter, we use nonrclativistic formulas having 
in mind that the order of magnitude of hadronic masses 
is m ~ 1 GeV while the freeze-out temperature T is of 
the order of 100 MeV. Note that even at the top SPS 
energy, where the pion yield almost twice exceeds the 
yield of all other particles, the system at the freeze-out 
stage consists mainly of heavy baryon and meson reso- 
nances which only later decay into pions. At this stage, 
thermal pions comprise less than 50% of all observed 
pions. 

From the kinetic theory of simple gases one can esti- 
mate the shear viscosity coefficient as [5lj 



havior displayed in Fig. [7J Here the experimental scal- 
ing turns out to be well reproduced in its high-density 
part, while the calculated elliptic flow still reveals no 
scaling at low charge densities. Note, however, that the 
lowest incident energy, at which the pion elliptic flow 
have been measured, is E\ & \, = 11.5^4 GeV (the E877 
data [3] in Fig. [7]) . If we consider vi only at high en- 
ergies £i a b > 10A GeV (lower bold lines in Fig. [7]), the 
reproduction of the scaling behavior becomes much bet- 
ter. In principle, if we let the fitting parameter c s ot r 
vary with the incident energy (i.e. smoothly join the 
corrected V2 at high energies with those kept unchanged 
at low energies), then the scaling at E\ a b > 10^4 GeV 
can be better reproduced. However, we deliberately 
avoid such a multi-parameter fit in order to keep the 
underlying physics more transparent. 

At lower AGS energies the 3FD model certainly pre- 
dicts no scaling. Apparently, this is a consequence of 
the partial shadowing of the transverse expansion by 
spectators. At lowest AGS energies this shadowing be- 
comes dominant and leads to the negative elliptic flow. 

It should be mentioned that scaling factors taking 
into account fluctuations of the initial eccentricity lead 
to a better quality of the experimental scaling [47| (at 
least at RHIC energies). Since the AGS-SPS data ana- 
lyzed here were obtained without applying the fluctua- 
tion corrections, we also do not use them in the scaling 
factors. Another aspect is influence of initial-state fluc- 
tuations on the hydrodynamic results. In our calcula- 
tion we do not introduce such fluctuations. In Ref. [48[ 
it was shown that the initial-state fluctuations can no- 



r)~-(p)n\. 



(6) 



where (p) ~ y^mr/V is an average thermal momen- 
tum of a typical hadron, A = (nctr) -1 is its transport 
mean free path and n is the total particle density of the 
medium 5 . For the rough estimate we neglect the angu- 
lar anisotropy of hadron-hadron cross sections. Taking 
otr = 6mb, obtained from the fit of vi data in Sect. Ill, 
we get the estimate 



(1.4-1.7) fm" 



(7) 



for T = (0.1 — 0.15) GeV. One can estimate the role 
of thermal pions (not hidden in resonances) by using 
the relativistic generalization of the above expressions 
suggested in Ref. [53[ . The direct calculation shows that 
pions contribute to 77 not more than 30% in the same 
domain of T. 

The relative strength of dissipative effects in fluid dy- 
namics can be characterized by the ratio of 77/s where 
s is the entropy density. The latter quantity can be 
estimated by using the Sackur-Tetrode formula [54[ 



In 



mT\ 3/2 1 



5 

+ -n. 



(8) 



5 The calculation of Ref. |52|1 for the hard-sphere gas gives a 
similar expression for r\ with <rt r = 47rrg , where ro is the radius 
of the sphere. 
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For T = (0.1 - 0.15) GcV and n = 0.2 fm~ 3 this for- 
mula leads to the estimate s ~ (0.6 — 0.7) fm~ 3 . The 
same calculation at n = 0.4 fm~ 3 gives higher values 
s ~ (0.8 — 1.1) fm -3 . However, Eq. ([5]) corresponds 
to the gas of identical particles and, therefore, under- 
estimates the entropy. The calculations for the gas 
of hadronic resonances with excluded volume correc- 
tions [551 ] gives entropy densities which are by about a 
factor of two higher. Finally, in the considered domain 
of T and n we get the estimate 77/s ~ 1 — 2 . At lower 
incident energies, noncentral rapidities and/or for more 
peripheral collisions the 77/s ratio may be even higher. 
This is expected, because the ideal hydrodynamics be- 
comes less justified in these regions. 

The obtained value of 77/s is of the order of mini- 
mal values observed in water and liquid nitrogen [56[. 
Therefore, the strongly-interacting matter in the con- 
sidered energy range indeed behaves like a viscous liquid 
rather than a gas. However, the obtained large values 
of the Knudsen number Kn ~ 0.7 may be interpreted as 
if the fluid dynamics is not applicable at all. This would 
be indeed so, if the estimated viscosity corresponded to 
entire stage of the hydrodynamic expansion. In fact, 
the estimated Knudsen number and 77/s ratio should be 
also attributed to late stages of the expansion after the 
freeze-out, i.e. to the afterburner. During this late stage 
the hydrodynamics becomes inapplicable and hence it 
is not surprising that the resulting effective Knudsen 
number turns out to be so large. This implies that 
Knudsen numbers corresponding to the hydrodynamic 
expansion may be still small, making hydrodynamics 
applicable. The argument in favor of such an inter- 
pretation is that the 3FD model is able to reproduce a 
large body of experimental data in the AGS-SPS energy 
range [IK . Another evidence in favor of this interpreta- 
tion follows from the study [1, H3] of a post-freeze-out 
evolution which has shown that the afterburner is re- 
sponsible for the major part of the discrepancy between 
the ideal hydro predictions and experimental data for 
i>2- Therefore, the above estimated values of 77/s should 
be considered as an upper limit for of the hydrodynamic 
expansion. 

The above arguments are not only relevant to the 
AGS-SPS energy domain. A similar analysis of the 
PHOBOS data at RHIC energies also gives rather large 
values of the Knudsen number. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have analyzed the elliptic flow in 
the SIS- AGS-SPS energy range within the 3FD model. 
Direct hydrodynamical calculations of the elliptic flow 
result in a good reproduction of experimental data at 
SIS and lower AGS energies (description of the FOPI 
pion data is still far from being perfect) but consider- 
ably overestimate the data at top AGS and SPS ener- 
gies. The latter problem cannot be cured by neither 



variation of the EoS stiffness nor the freeze-out crite- 
rion. Changing the inter-fluid friction forces does not 
solve this problem either. In this paper we attribute this 
problem to dissipative effects during the expansion and 
freeze-out stages. In order to estimate the role of dissi- 
pation from the difference between the 3FD results and 
the observed data, we apply a simple formula, proposed 
in Ref . (30j , where this difference is expressed in terms 
of the Knudsen number. It is shown that the interpre- 
tation of the disagreement between the 3FD model and 
experimental vi data due to dissipation turns out to be 
fruitful. With the help of a singe parameter we are able 
to fit the calculated tj 2 values to the observed data in 
a broad range of incident energies, ccntralitics and ra- 
pidities. Moreover, the experimental scaling behavior 
of pion elliptic flow scaled with the initial eccentric- 
ity versus chargcd-hadron-multiplicity density per unit 
transverse area turns out to be reasonably reproduced. 

Proceeding from values of the Knudsen number, de- 
duced from the above fit, we estimated the viscosity-to- 
entropy ratio as 77/s ~ 1 — 2. These values of 77/s are of 
the order of minimal values observed in water and liq- 
uid nitrogen [56| . Therefore, the nuclear matter in the 
considered energy range indeed behaves like a liquid, 
however, not so perfect as at RHIC energies. 

The estimated 77/s ratio accumulates the effects of 
dissipation at the hydrodynamic expansion stage and 
the afterburner stage after the hydrodynamic freeze- 
out. Having in mind that the afterburner may give an 
important contribution to dissipation we conclude that 
the estimated 77/s values represent only an upper limit 
for this quantity when applied to the hydrodynamic ex- 
pansion stage until the kinetic freeze-out. On the other 
hand, the authors of Ref. [53| have estimated the 77/s ra- 
tios on the chemical freeze-out line observed for central 
collisions of heavy nuclei at beam energies from SIS to 
RHIC. Within the excluded- volume hadron-resonance- 
gas model, they found that 77/s are larger than 0.3-0.5 
depending on the hard core radius. We would like to 
stress that these values should be lower than ours be- 
cause they correspond to the earlier chemical freeze-out 
stage. 

Recently the CERES collaboration reported that 
77/s = 0.021 ±0.068 at the top SPS energy [57|. This re- 
sult is certainly in contradiction to our estimate of this 
ratio. In Ref. [57J the 77/s ratio was deduced from anal- 
ysis of two-pion correlation measurements, more pre- 
cisely, of the transverse-momentum dependence of the 
longitudinal pion source radius. The analysis was per- 
formed on the basis of boost-invariant (Bjorkcn) solu- 
tion to the relativistic viscous hydrodynamics. Con- 
cerning the latter we would like to mention that it is 
not well justified experimentally since the pion rapidity 
spectra are far from being flat at SPS energies. As it is 
illustrated in Ref. [58[ , the 3-fluid dynamics of nuclear 
system at 158^4 GcV is quite different from the Bjorkcn 
picture. This could be a reason of disagreement between 
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our results and the results of Ref . [57} ■ 
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